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Abstract 

We demonstrate by mathematical analysis and systematic computer simulations that redistribution can 
lead to sustainable growth in a society. In accordance with economic models of risky human capital, 
we assume that dynamics of human capital is modeled as a multiplicative stochastic process which, in 
the long run, leads to the destruction of individual human capital. When agents are linked by fully- 
redistributive taxation the situation might turn to individual growth in the long run. We consider that 
a government collects a proportion of income and reduces it by a fraction as costs for administration 
(efficiency losses). The remaining public good is equally redistributed to all agents. Sustainable growth 
is induced by redistribution despite the losses from the random growth process and despite administrative 
costs. Growth results from a portfolio effect. The findings arc verified for three different tax schemes: 
proportional tax, taking proportional more from the rich, and proportionally more from the poor. We 
discuss which of these tax schemes performs better with respect to maximize growth under a fixed rate 
of administrative costs, and with respect to maximize the governmental income. This leads us to some 
general conclusions about governmental decisions, the relation to public good games with free-riding, and 
the function of taxation in a risk taking society. 

Introduction 

This paper shows how redistribution of income spurs growth of human capital in a society just because 
of a portfolio-effect. Our model captures the concept of "risky human capital" from recent economic 
literature (see [l]), where human capital is described by a multiplicative stochastic process which, in 
the long run, leads to the destruction of individual human capital. We model the random growth or 
decline of human capital as proportional to the individual endowment of income. We couple agents by 
fully-redistributive taxation [2] (that means collected taxes are equally redistributed) which is associated 
with efficiency losses j3j. Redistribution re-balances gains and losses from individuals and works as a 
portfolio effect which spurs growth into the individually lossy stochastic processes. 

Economic literature does not explicitly point out the portfolio effect through redistribution in the 
relationship between inequality and growth of human capital. So far, in the politico-economic literature, 
three basic reasons are discussed of why redistribution is beneficial for society. The first branch of 
literature stresses the insurance aspect through redistribution (see [i]). Mirrlees highlights that from a 
welfare maximizing point of view, the level of redistribution should be at a level on which the poor do not 
suffer and both the poor and the rich have an incentive to improve their situation |4]j^] More recently, 
the literature of socio-political unrest (7lp| points out that redistribution guarantees social stability and 
reduces the effort the society has to make when inequality is high. 

1 Present address: Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany 

2 Mirrlees motivates the insurance aspect with fairness-considerations. Therefore, this literature is also partly linked to 
the literature about social preferences [5][6]. In a series of experiments, it is shown that subjects have a preference for 
avoiding high degrees of inequality even if they have to resign payoff. 
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In the second branch, redistribution reduces the disincentive for the poor for taking too high risks [9] . 
Redistribution increases the endowment of the poor. The poor reduce the demand for loans and invest 
more efficient which means that they take less risk. The efficiency in the economy is improved, growth 
will be higher. Redistribution spurs growth because of decreased disincentives in society. 

The third branch describes the transmission channel of the median- voter-approach between inequality 
and growth [3]. The utility maximizing calculus of the median- voter determines the level of redistribution. 
The median voter proposes the median level of redistribution. This level is unbeatable in pairwise majority 
decisions. If the inequality is high, the median- voter enforces a high level of redistribution. Redistribution 
is assumed to induce efficiency losses, as in our model. High levels of inequality are associated with high 
efficiency losses and therefore low growth rates. For instance in [lOj it is shown that high degrees of 
inequality are not associated with high levels of social spending. This motivated the modification and 
extension of this approach by the consideration of institutions [If] and elites |12| . 

In contrast to above mentioned work on the relationship between inequality and growth, we argue that 
neither socio-political nor incentive considerations have to be taken into account to show that redistri- 
bution spurs growth. Our model excludes incentive-incompatibilities, voting-approaches and normative 
concepts like insurance or fairness considerations, to point out the effectiveness of the pure portfolio effect 
through redistribution. 

As we will see, the effect of combining lossy proportional stochastic growth and linear lossy redistri- 
bution is non-trivial as both processes lead to destruction of human capital when they run independently, 
while their combination can enable survival. This "magic" effect of induction of growth from two lossy 



processes is based on the portfolio effect known from investment science 13 : Gains and losses are re- 
balanced by redistributing income into human capital "assets" , which ensures optimal growth of the 
portfolio. The effect has been discussed before under different names, such as repeated Kelly games, 
Kelly optimal portfolio, and re-balancing of asset allocations, and seems to be rediscovered from time 
to time in a new context 14-22 . It's applicability in various fields is laid out in 123]. We will call the 
phenomenon portfolio re-balancing effect in the following. This interprets each subject with its human 
capital endowment as an asset and the society as the portfolio. Asset values change stochastically, and 
taxation and redistribution of income re-balance the values of the different assets. 

The econophysics literature has studied extensively the statistical mechanics of money based on its 
conservation 24 ■ 26 and within a kinetic exchange model 27 (see also the paper and the literature review 



in [28]). The topic of taxation and redistribution has been introduced in these models [29]. Regarding 
this stream of literature, we focus on the growth enhancing effect of redistribution neglecting conservation 
of money. 

The purpose of this paper is to discuss the role of the portfolio re-balancing effect of fully-redistributive 
taxation for the growth of society's human capital. By systematic computer simulation we quantify 
conditions of the stochastic growth process, and the taxation schemes that prevent the destruction of the 
accumulated human capital of the society. Further on, we quantify optimal tax rates which maximize 
growth of society's human capital, and how a selfish government optimizes its income. 



Methods 

We describe our method by the following steps: 

• Specification of the outlined economic model on the endogenous variables human capital h and 
income y in terms of a straight forward transition of human capital to income, redistribution 
of income between agents and independent random multiplicative production of human capital 
proportional to income. Exogenous variables are: tax rate a, rate of administrative cost 6, taxation 
scheme, the distribution of random growth factors, and the number of agents N. 

• Theoretical demonstration what conditions are most interesting for studying the effect of portfolio 
re-balancing on growth. 
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• Theoretical demonstration of dynamics at border cases. 

• Presentation of some example runs of the process. 

• Description of the setup of a systematic simulation the extraction of the average growth factor from 
simulation data. 



Specification of the economic model 

Let us consider a society of N agents, each of which is characterized at time t by its human capital hi(t) 
which is a positive scalar value. Consequently, h(t) € M. N is the human capital vector of all hi(t). Let us 
denote the total human capital at time t by H{t) = 5^ i=1 hi(t). 

The production process is: human capital is used to produce income yi income is taxed and fully- 
redistributed, and income is directly invested in human capital. This is formalized by the equations: 

y i (t) = W od i (h i (t)), (1) 
W (t+l) = rediSi(tf(t)), (2) 
fti(t +l) = HCprod i ( W (t+l)). (3) 

Notice that prod^ and HCprod i are functions operating on individual values, while redis^ is scalar-valued 
but takes the whole income vector as input. These scalar- valued functions serve as component functions 
for the vector- valued functions prod, HCprod, and redis which are self-maps on l w . 

Given an initial human capital vector h(0) the evolution of human capital in the system of individuals 
is described by the equation 

h(t + 1) = HCprod(redis(prod(/i(i)))). (4) 
Consequently, the evolution of income is given by 

y(t + 1) = reois(prod(HCprod(y(i)))). (5) 

We call Y(t) — J2iLi Vi(f) tne total income. Growth (positive or negative) of this aggregated variable is 
analyzed in the following. The total income evolves of course tightly related to the evolution of the total 
human capital H(t) because of Eq. 



Production of income and human capital We assume that production is directly transfered into 
income 

y L = prod^) = hi, (6) 

which means that the wage is equal to one. The production of human capital is assumed to be based on 
an individual multiplicative stochastic event 

hi = HCprod, (j/i) = rn{t)yi (7) 

where rji(t) is a realization of the positive random variable r\. If agent i at time t has income yi(t) 
then after producing its human capital is rji(t)yi(t). When rji(t) < 1 human capital declines, otherwise 
it grows. Thus, we assume that human capital for the next round of production is build from current 
income times a random factor. This can be interpreted as a generation model where each generation lives 
for one period and invests its income in the human capital of its successor. In an innovative economy 
the value of old human capital quickly dissolves and new human capital must constantly be produced by 
investing income. Thus, the human capital production function is also reasonable on shorter time scales 
then generations. 

Without subscript i, prod and n(t) are meant as vectors. Thus, growth dynamics of human capital 
vectors h reads HCprod(?/) = r](t)y. The product T)(t)y is meant component- wise, r)(t) being an equally 
sized vector of independent realizations of n. 
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Our production function has only the input factor human capital and is therefore simplified. The 
simplification is motivated by the idea of the human capital intensive production in modern economies 
which is in line with endogenous growth theory (see |30| ) . For reasons of comparability, we present in 
the following some feasible extensions of our model. A standard Cobb-Douglas production function also 
includes the input of capital k and labour I (accompanied by their exponents < a, f3 < 1 with a + /3 = 1). 
Moreover, the transfer of income into human capital is accompanied by consumption c and saving rates 
s. This would imply Equation ^ to look as 

Hi = prod(/ij) = hik a l (i 

and Equation ^ to be 

hi = HCprod( yi ) = tfc(i)(l - s)(l - c) Vl . 

Equation ^ would read 

y(t+l) =redis( m {t){l - s){l - c) Vl k a f ). (8) 

When s, c, k, I, a, (3 are all constants the term rji(t)(l — s)(l — c)fc a Z' 3 is a draw from a random variable 
from the distribution as rj but scaled by a constant factor (1 — s)(l — c)k a l@ . Modification of these factors 
have thus the same effect as a multiplicative scale in the random variable rj. By holding k constant, we 
implicitly assume that savings are equal to depreciation of capital. Our model is thus based on a standard 
economic growth model and is simplified to focus on the portfolio re-balancing effect. 



Redistribution We quantify the redistribution function with three different taxation schemes: propor- 
tional taxation, a progressive scheme where agents have to pay everything above a dynamically chosen 
maximal tax-free income, and a regressive scheme where agents have to pay either a dynamically chosen 
fee - like a per capita premium - or all their income if they cannot afford the full fee. (No worries, 
agents get back some income because of the redistribution.) All three schemes we specify by the same 
two independent parameters: the tax rate a, which determines the fraction withdrawn from the total 
income of all agents, and the rate of administrative cost b, which determines the fraction withdrawn by 
the government from the raised taxes before redistribution to agents. 

Let us call the amount of taxes collected from agent i to be taxj(y). Notice, that it depends on 
the vector of income. This enables us to define dynamically adjusted taxation schemes which take the 
distribution of income into account. Naturally, taxj(y) should be confined between zero (no tax) and y^ 
(tax equals income). 

The tax revenue is collected by a government at a central place, which involves administrative costs 
(efficiency losses cf. [3]), which are assumed to be proportional to the amount of taxes raised, i.e. b e [0, 1] 
denotes the rate of administrative cost. Consequently, the public good for redistribution is the raised taxes 
minus the cost: 

N 

pg(y) = (l-6)^tax i (y), (9) 
i=i 

while the government income is 

N 

gi(y) = 6^tax l (y). (10) 

»=i 

Because of a fully-redistributive tax, pg(y) is divided with equal shares among all agents, i.e. for every 
agent income increases by an amount pg(y)/N. Other mechanisms of redistribution in a related model 
are analyzed in [29] . The redistribution function (net income) for agent i is thus 

rediSi(y) = y t - taxi(y) + ^jr~- ( n ) 
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Depending on the position within society, an agent could be a net tax payer or a net transfer recipient. 
We specify the tax function taXj(y) with respect to the tax rate a £ [0, 1] such that it holds 

^taxi(y) = a^yi = aY. (12) 

i i 

We distinguish three schemes of taxation which differ in from whom the fraction a of the total income 
is raised: (i) proportionally from everyone, (ii) more than proportionally from the poor (regressive), or 
(iii) more than proportionally from the rich (progressive). In progressive and regressive taxation schemes 
tax rates differ in given income brackets. We consider extremal cases to pronounce differences. 

(i) Proportional tax is the classical taxation scheme where each agent has to pay a fraction a of its 
individual income 

taxj(y) = ay { . (13) 

(ii) Regressive tax charges a fixed fee Cf cc (y) > from everyone if possible, otherwise all income is 
charged. 

taxi(y) = min{x i,c icc (y)} (14) 
In the latter case, the agent still receives its proportion from the public good, so it will not be 



without income after redistribution. The fee has to be such that ( 12 ) is fulfilled for the current 
income vector. (It is easy to see that this is possible and unique.) 
(iii) Progressive tax charges all income exceeding a threshold of tax free income c max (?/) > from every- 
one. Every agent with income below c max (?/) > pays no taxes. 

taxj(y) = maxly, - c max (y),Q} (15) 



The threshold has to be determined such that ( |T2| ) is met. (It is easy to see that this is possible 
and unique.) Notice, that the income of an agent who has to pay taxes is larger than c max (y) after 
redistribution because of its share from the public good. 
The three schemes are comparable in that the total amount of raised taxes is always a fraction a of 
the total income, regardless of the shape of the distribution of income. Thus, they all deliver a public 



good of pg(y) = (1 — b)aY . Note, that in the related model of 29 taxes are raised when agents trade 
and not every time period from every one. 

The regressive and the progressive tax assign different tax rates in two different tax brackets [0, c(y)] 
and [c(y), oo]. The regressive tax scheme (where c = Cf co ) taxes 100% in the lower tax bracket and 0% in 
the upper. The progressive tax scheme (where c = c max ) taxes 100% in the upper tax bracket and 0% in 
the lower. Figure [I] demonstrates an example with six agents of different income. It is shown what will 
be charged from each agent and how income looks after redistribution for each of the three schemes. 

Note, that the dynamic fee Cf CB (y) and the dynamic maximum c max (y) are implicitly defined, to meet 



the condition of Eq. ( 12 ). In realistic taxation systems, it might seem impractical to determine the fee and 
the maximum after the current income of all agents is known. In reality, one would only adjust thresholds 
for the next turn. We omitted that modification to prevent delay effects. Probably, this modification 
would cause only minor changes. 



On the distribution of random human capital growth factors 

Analysis of human capital production without redistribution (a = in Eq. (j2J)) does not involve interac- 
tion. Therefore, it is enough to focus on a single agent and Eq. ^ collapses to y(t + 1) = rj(t)y(t). Let 
■q have finite variance. With y(0) = 1 it holds 

y(t + 1) - V (t)y(t) = v(t)v(t -!)••• v(X)v{0) = J[ 17(a). (16) 

s=0 



() 



This resembles the human capital life cycle model [T] where human capital h t at time t is determined as 
the result of a stochastic process h t = (Tt-i ■ ■ ■ o\Qi where i is the initial investment in human capital, 9 
is a stochastic shock to the human capital investment and a s are stochastic human capital depreciation 
shocks. 



Eq. ( 16 1 is equivalent to 

t 

log y(t + 1) = log V {t) + log y(t) =J2 lo § V(s). (17) 



The central limit theorem applied to (17) implies that the distribution of the random variable logy(t) 



gets closer and closer to a normal distribution J\f((i t ,a t ) with mean and variance parameters 

IH = *Mlog7) 5 <?\ = Wogr, ( 18 ) 

with /iio gJ) = (log rj) and &\ ogrj = ((log?/) 2 ) — (log?/) 2 . Consequently, for t — > oo the distribution of 
y(t) approaches the log-normal distribution \og-J\f(t/Ji ogri ,\/tai g rj ). Based on that fact, we chose the 
log-normal distribution with its two characterizing parameters as the distribution of rj in our simulation 
setup. 

The expected value of income might grow, while every individual trajectory of y(t) dies out. The 
condition for this seemingly contradictory situation is 

Hiogt) < < log/j r) (19) 

which is equivalent to (??) goo = exp/ii ogt? < 1 < /i^ = {rj) i.e. the arithmetic mean of r\ being larger 
than one, while its geometric mean is less than one. Elementary explanations of this effect are given 



in 13 22 31 . It can be shown that for long enough time span any single trajectory grows only with the 
geometric mean (??) goo . For log-normal distributions of rj the two inequalities in (19 1 are equivalent to 



Wogr; being negative while crf ogI) is sufficiently large af ogv > -2/i logJ? . 

This situation forms the basis for the effect of growth which is induced by coupling lossy multiplicative 
stochastic growth with lossy redistribution in a finite population. Redistribution helps the system to 
realize a growth rate somewhere in between the geometric and the arithmetic mean of rj. 

To reduce the number of independent parameters in simulation, we choose log-normal distributions 
where (r/) ■ (r)) aco = 1 holds. Under this condition, the two parameters of the log- normal distribution \x 
and a (mean and standard deviation of the underlying normal distribution) are represented by one free 
parameter which allows for different skewness, but keeps the balance of the expected value (rj) an the 
realized growth rate (??) goo - This condition enables that destruction and growth are theoretically possible 
for distributions of this class. 



Theoretical analysis of border cases and an example 

We are interested in the average realized growth rate and its dependence on the independent parameters. 
For some border cases we can theoretically derive that average growth is exponential Y(t + 1) = gY(t) 
and also quantify the magnitude of the growth rate g. 

Case 1: Only redistribution, no stochastic production of human capital ((?7) gco = 1 = (??))• For any 
taxation scheme, and any N it holds g = (1 — ah). Thus, there is never growth of human capital. 

Case 2: Only stochastic production of human capital but no redistribution (a — 0). For any taxation 
scheme, any b and any TV it holds g = (r)) Keo - 

Case 3: 100% tax and infinite number of agents (a — 1, N = oo). All trajectories act as one, and for 
any taxation scheme g = (1 — ab) (r/). The growth with the mean can be realized because in an infinite 
society no rare but large growth event is "missing" . 
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Based on the last two cases we argue, that for intermediate a and finite N the average growth rate 
lies somewhere in between, but we do not have an analytic expression for it. 

Figure |2j gives an example, where we fix (77) = 3/2 and (i]) gco — 2/3. This implies log-normal 
parameters /1 = —0.405 and a = 1.274, which we use as the distribution of intermediate risk in simulation. 
Under this distribution income declines with a probability of 62.5%, it at least doubles with probability 
19.4%, and it will be more than ten times larger with probability 1.7%. Tax and admin rates are set at 
intermediate levels a = 0.3, b — 0.2. Trajectories are computed according to ^ for a society of N = 10 
agents, each starting with human capital equal to one. Trajectories are shown for all three tax schemes. 
Each trajectory is computed with the same random draws from the random variables r]i(t) for each i and 
t, to allow for a direct comparison of the different taxation schemes. 

In Figure [2] progressive taxation (iii) leads to the largest growth of total income, proportional taxation 
(i) resulted also in a growing society while regressive taxation (ii) fluctuates between growth and decline 
with no clear trend visible. When there was no redistribution at all there would be decline (with a growth 
factor (r]) gco = 0.667). Pure redistribution without production of human capital would also imply decline 
(with a growth factor of 1 — ab — 0.94). The performance ranking progressive better than proportional 
better than regressive holds even if we vary the essential parameters admin rate 6, tax rate a and the size 
of the population iV as shown in the lower part of Figure [2] Regarding their impact on the dynamical 
behavior, we see that for a higher admin rate (b = 0.6) growth might turn to decline. A lower tax rate 
(a = 0.01), can also imply destruction of income and human capital for all three taxation schemes. In 
this case the portfolio re-balancing effect of redistribution is not used well and this is not compensated 
by the savings from the lower loss of redistribution (1 — ab = 0.998). Finally, in a larger society with 
N = 100 all schemes achieve larger growth factors. 

Progressive taxation may contribute to disincentives for the decision to invest income into the pro- 
duction of human capital. Especially in our extremal case income above a maximal income is taxed by 
100%. To that end, let us assume that every agent can decide what fraction of its income to invest while 
the remaining income's value remains as is. Taxation and redistribution remains obligatory. If the agent 
is to maximize the value of its income after human capital production, production and redistribution, 
the rational decision would be to invest all, as long as the expected value of human capital is larger 
than the invested capital. This holds also for progressive taxation, as any increase of expected value 
before redistribution increases the expected value of the public good and thus the own expected value 
after redistribution. In our model, progressive taxation does not remove the rationality of the choice of 
investing all into production of human capital, but makes difference to not investing in human capital 
smaller. We do not touch the issue of free-riding (which would be to avoid paying taxes) in the above 
argument. When free-riding was possible, not paying taxes but receiving a share from the public good, 
is of course rational under every taxation scheme. 

Simulation setup, independent and dependent variables 

For each of the three taxation schemes we aim to get an overview about the dependence of the average 
growth factor on the tax rate and the rate of administrative costs. Further on, we want to control for the 
effect of society's size and lognormal distributions which are more or less risky (in the sense of higher or 
lower right-skewness). 

Consequently, we set up a systematic computer simulation to estimate the average growth factors g. 
Table [T] shows the list of independent variables, their ranges in the simulation setup. Growth factors are 
estimated on the basis of stochastic trajectories of total income Y(t) for 500 time steps by regressing 
log(g) in \ogY(t) — log N + t ■ log(g). The intercept in the regression is naturally fixed at log N. From 
100 estimated growth rates of such stochastic trajectories we compute the average growth factor as 
the geometric mean. The geometric mean is used because it better fits the central tendency of the 
distribution, as a growth rate is naturally a parameter larger than zero an thus typically log-normally 
distributed. Changing to the arithmetic mean would shift the results a bit, as the arithmetic mean is 



8 



always larger than the geometric mean. 

To get on overview about the parameter space we cover the (6, a)-plane by a fine grid while the 
number of different sizes of the society and different distributions was kept low to make computations 
finish within less than two weeks on a laptop. See the matlab-code in the supporting material which 
produces the simulation data (see function dataMSPgrowthrates). 

As our focus is on maximizing growth factors, optimal tax rates, and government income let us define 
further variables which depend on the average growth factor as a function of a and b, g(b,a): For a 
fixed admin rate b we define the maximal growth factor g max (b) = max a g(b, a) and the optimal tax rate 
Oopt(fr) = argmax a g(6, a). (Note, that argmax is not necessarily unique, but empirical results support 
the conjecture that there is only one local maximum, and consequently a op t(&) is unique.) The rate of 
government income at the current total income is gov(6, a) = baY(t + 1)/Y(t) — bag(b,a). The rate of 
government income under optimal tax rate as a function of b is defined by gov opttax (6) = ba op t(b) g ma ,x(b). 

Results 

Description of figures 

Figures [l}|4] are the core Figures to understand the message of our paper. Figures [l] and [2] explain and 
demonstrate our economic model and Figures [3] and [4] show main results. 

We visualize our simulation results exemplarily for N — 10 and the distribution of intermediate risk 
((n) , (?/) gco ) = (1.5,0.667) in Figure [3j We show in Fi gure 3A the average growth factor g(b,a) color- 
coded in the (b, a) parameter plane for each of the three tax schemes. In each plot, the solid line divides 
the "zone of sustainable growth of income" (yellow to red) from the "zone of income destruction" (yellow 
to blue). The dashed line shows the optimal growth maximizing tax rate a op t(b) for a given admin rate. 
In Figure [3j3, we show the critical lines dividing the zones of growth and destruction and the optimal 
tax rates in one plot to compare the three taxation schemes. (In all plots the black dotted line shows the 
maximally possible size of the zone of growth, where (1 — ba) (rf) = 1. Above this line it is trivial that 
income can not grow. 

Fi gure |4]A shows the maximal growth factor j ma x(^), |4]A the optimal tax rate a pt (6), and||C the 
rate of government income at optimal tax rate gov opttax (6). All are functions of the admin rate b and 
they are shown for all three taxation schemes in our standard color-code. 

Figure [5}|8] are extensions of Figures J3] and HI by showing also the data for N = 100, a less risky, and a 
more risky distribution of r\. Figures [5j |6j and[8j are the pure extensions of Figures [3]A.,[3j3, and|4j while 
Figure [7] is a regrouping of lines to different subplots focusing on comparison of JV = 10 and N = 100. 

On growth and destruction of income and human capital 

Figure [3] and related Figures [5j [6j and [7] summarize what combination of tax rates and admin rate allow 
for society's income and human capital to grow sustainable. It is interesting to focus on situations where 
cither the tax rate a or the admin rate b is fixed: 

Constant tax rate a: Raising the admin rate b always lowers the growth factor g{a,b) and turns 
the growth regime at some point into the destruction regime. 

Constant admin rate b: The average growth factor g(a, b) is not monotonic in a. For very high 
and very low tax rates, the growth factor is the lowest and can lead to income destruction, while only 
intermediate tax rates prevent this. High tax rates tend to lower the growth factor because a larger 
fraction of the total income is reduced by the admin rate (see the definition of the public good ([9])). On 
the other hand, very low tax rates lower the growth factor because the portfolio re-balancing effect is not 
used optimal, thus part of the human capital "gambles it self away" . 
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This characterization is ubiquitous for all taxation schemes, both population sizes and the riskier and 
less risky distribution of r/. It holds ubiquitous that the zone of sustainable growth of regressive taxation 
is contained in the zone of growth of proportional taxation, which is further contained in the zone of 
growth of progressive taxation, when N and the distribution of r\ are kept constant. When the taxation 
scheme and the distribution of r\ is kept constant the zone of growth of the smaller society (N = 10) is 
always contained in the zone of growth of the larger society (N = 100). 

These findings suggest that the portfolio re-balancing effect is used more effectively, when the society 
is large and when proportionally more is taken from the rich (progressive tax), than from the poor 
(regressive tax) . Simple inclusions of zones of growth do not hold for comparisons of low risk, intermediate 
and risky distribution of n. We refrained from comparing them in detail, because our balancing condition 
(77) • (r)) gco = 1 is somehow arbitrary. 

In the following we answer five questions about optimal choices of tax and admin rates for different 
perspectives. 

On growth maximizing tax rates and taxation schemes 

(a) What is the optimal tax rate a op t(b) and how does it differ between the three taxation schemes? The 
optimal tax rate is 100% under admin rate b — for any tax system but it declines fast with rising admin 
rate as can be seen in Figure |4]A.. Within the range < b < 0.35 the progressive tax scheme reaches 
the lowest optimal tax rate, regressive taxation the highest. For a realistic admin rate of about 20% the 
optimal tax rate in the progressive tax scheme and the proportional tax scheme is less than 30%, but it 
is larger than 50% under the regressive tax scheme. The ranking is inverted for large admin rates. From 
Figure [8] it can be seen that these rankings and the switch of the ranking also holds for larger societies 
being more drastic for low admin rates and less drastic for high admin rates. For riskier societies the 
optimal tax rates are larger in general. 

(b) Which taxation scheme reaches the largest average growth factor for a given admin rate and 
optimal choice of the tax rate? As can be seen in the central panel of Fig.|4jB the progressive tax scheme 
achieves the highest maximal growth factors for all admin rates. The proportional tax scheme is always 
second and the regressive tax scheme ranks last. Hence, the largest growth factor is reached with the 
progressive tax scheme that takes more than proportionally from the rich. 

On income maximizing governments 

Let us assume, that governments are forced to choose tax rates close to the optimal tax rate. A government 
might be forced to do so in an informed and democratic society, when the impact on the tax rate on 
average growth is known and when voters wish that growth rates of total income or human capital are as 
large as possible. Under this assumption the rate of government income under optimal tax rate gov opttax 
is of interest, because we can ask what admin rate a government might choose to maximize its income. 
The rate of administrative costs is usually also under the control of the government, but we assume that a 
government is not forced to optimize it for growth. Optimal would of course be to have no administrative 
cost. One reason is that it might be easier for a government to argue, that the administrative cost 
cannot be lowered, because of fixed contracts. Another argument is that an alternative party which could 
overtake government has possibly the same interest of increasing its income and consequently democratic 
competition does not work as easy. Based on these two assumption and simulation results we can answer 
three further questions: 

(c) Which admin rate b would a self interested government choose? Figure [4]C shows that the rate 
of government income under optimal tax rates is not monotonic in b. In particular, for higher values of 
b the growth of total income becomes smaller, hence even a government maximizing its income has no 
incentives to raise the admin rate to the largest possible. This is because large admin rates reduce growth. 
The admin rate where the government income is maximal is marked by "*" in all three panels. It varies 
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with the taxation scheme: lowest admin rates for regressive tax, highest admin rates for progressive tax. 
This ranking only changes for the riskier distribution of r\ (see Figure [8]) , with proportional tax having 
lowest admin rates. 

(d) Which taxation scheme would a self interested government choose ? Looking at absolute values of 
gov opttax (&) a self interested government would choose the regressive tax because it gives the maximum 
income of all schemes, even at moderate admin rates. Thus, the largest government income is reached 
with a scheme that takes more than proportional from the poor. This is caused mainly due to the fact 
that the optimal tax rate is much larger under the regressive tax scheme. Consequently, the share raised 
by the admin rate is also larger as under other schemes. Thus, this result crucially depends on the 
assumption that a government is forced to implement the optimal tax rate, while the taxation scheme is 
given. 

(e) Which taxation scheme delivers the largest average growth factors under a self interested govern- 
ment? Looking at the "^"-symbols in Figure |4p which come from optimal admin rates with respect to 



gov opttax (6) in Figure |4]C, we find that proportional tax delivers the highest average growth factors. This 

holds although regressive taxation can deliver the highest income for the governemt as seen in (d). The 
reason is that the regressive tax has much lower growth rates than the other schemes in general. It also 
holds although progressive taxation always delivers higher growth rates than the other regimes as seen 
in (b). The reason here is that the progressive tax attracts the government to raise the admin rate to 
optimize its income. 

Discussion 

In our view, redistribution enhances the dynamic potential for human capital production of an economy. 
The enhancement can be explained by the effectiveness of the portfolio effect. The enhancement seems 
is possible for proportional, regressive and progressive taxation schemes. 

The answers to questions (b) , (d) , and (e) deliver three different choices of one of the three taxation 
schemes, (b) suggest that the progressive taxation scheme should be chosen because it is always superior 
when a rate of administrative cost is fixed. Consequently, taking proportionally more from the rich is 
socially optimal when we can assume that the rate of administrative costs is an externally fixed parameter, 
(d) instead suggests, that a selfish government would decide for the progressive taxation scheme because 
under this scheme growth optimizing tax rates are much higher which leads to higher government income. 
Finally (d) shows that from the three taxation schemes proportional taxation achieves the highest average 
growth factor under an income maximizing government which can freely adjust the rate of administrative 
cost. The regressive scheme turns out to be too inefficient in turning on the portfolio re-balancing effect 
to enhance growth, while progressive taxation turns out to give incentives to raise the admin rate to such 
an extend that the loss due to this outweighs its efficiency in using the portfolio re-balancing effect. 

By focusing on the portfolio re-balancing effect we proposed a new approach to think about the link 
between inequality, redistributive taxation and growth. With our simulation, we have shown that taxation 
and redistribution can be a crucial ingredient to ensure the survival and development of a society relying 
on risky multiplicative stochastic growth of human capital. Our approach gives another explanation 
of why redistribution can be beneficial for growth. Together with the other approaches mentioned in 
the introduction, we show that the interplay between inequality, redistribution and growth depends on 
preferences (fairness and insurance considerations), socio-political unrest (social stability), incentives and 
disincentives (effort), the calculus of the median voter (voting-system) and on the portfolio re-balancing 



When paying taxes was voluntary, payment could be seen as an act of cooperation, which seems 
irrational but ensures the sustainable growth of human capital. As in the related public goods game, 
this society would be vulnerable to free-riders. In the public goods game the free rider problem is often 
solved by social norms or governmental forces to pay taxes. 




effect. 
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At difference to the classical public goods game, where the public good is multiplied by an efficiency 
factor larger than one, our model does not have such an amplification. On the contrary, from the collected 
public good a fraction is subtracted for administrative costs, which is equivalent to an efficiency factor 
less than one. Consequently, the emergence of cooperation, i.e. the sharing of income in order to sustain 
a long term growth, is even more subtle in our redistribution model, because even a normative call "Pay 
your part and it will be immediately increased by an efficiency factor!" does not work as easily. 

If we assume that different societies compete, evolution would promote those societies with higher 
overall growth factors of their total income. Thus, there should be an evolutionary adaptation towards 
the optimal tax systems without assuming other forces. Such an idea is closely related to group selection 
as a mechanism to promote the evolution of cooperative behavior [32] . Perhaps, the portfolio re-balancing 
effect is also a reason for the evolutionary success of religions which propose something like the tithe. 
Tithing 10% of income to charity in a religious society ensures better growth of society's human capital 
which might be an evolutionary advantage against other societies because of the portfolio re-balancing 



effect. The re-balancing effect might also be of relevance in other areas, such as biodiversity 33 or 
knowledge sharing, to enhance innovativeness in social and economic systems. 

Can we draw conclusions for large societies of some millions as in real world societies? We speculate 
that sizes larger than TV = 1000 imply even lower optimal tax rates because the portfolio re-balancing 
effect works even with very low tax rates and thus admin costs can be saved by low tax rates. But on 
the other hand we speculate that riskiness of individual stochastic growth also rises in larger societies 
which consequently implies higher optimal tax rates (see Figure [8]). It is much less likely to gain the best 
fitting human capital in a large society because there are more competitors. On the other hand, having 
gained the right human capital might lead to large benefits because there are many customers which 
could benefit from it. In conclusion, we speculate that our results are probably still valid for societies of 
real world sizes. 
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Table 1. Simulation setup. 



independent variable 


range 


size of range 


tax rate a 


+.°.-9 2 1 


51 


admin rate b 


0,+.°.? 2 ,0.8 


41 


taxation scheme 


regressive, proportional, progressive 


3 


riskiness 


less risky intermediate more risky 




((V) > ( J 7>geo) 


(1.25, 0.8), (1.5, 0.667), (3, 0.333) 


3 


number of agents N 


10, 100 


2 


^max 


500 


1 



Numbers of parameter values multiply to 37,638 combinations. 100 simulation runs with hi(0) = 1 were 
computed for each. Consequently, the growth factor g was estimated regressing log(g) in 
logY(t) = log A + t ■ log(g) (notice the intercept is naturally fixed at log A). For these 3,763,800 values 
of g the geometric mean was computed for each combination over all 100 runs giving 37,638 average 
growth factors as the basis for Figures [3}]8] 
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Progressive Tax 
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Figure 1. Demonstration of the redistribution function. Proportional tax, regressive tax and 

progressive tax functions. Six agents with income y = [100, 300, 600, 1000, 1500, 2100]. All tax functions 
are such that a — 1/3 of the total income (= 5600) is taxed, the administrative cost is set to 25% 
(b = 0.25). In numbers: pg(y) = 1400, which implies Cf cc = 366| and c max = 911 ^ . 
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N = 10, a = 0.3, b = 0.2 




200 400 200 400 200 



Figure 2. Example trajectories for tax schemes and exemplary modifications in subfigures 

Trajectories computed with the same realizations of the random variables rji(t) and iteration of Eq. ^ 
Black lines show limiting cases: lower line shows "no tax" Y(t) = N({<q) ) = 10 ■ (0.667)*; upper line 
shows "full tax and infinite number of agents" Y(t) = N((l - ab) (77)) = 10 • (1.41)*. 
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color axis: average growth factor g(b,a) 
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Figure 3. Zones of growth and destruction. Based on the average growth factor g(b, a) for 
N = 10, (r/) = 1.5, and (?/) gco = 0.667 (cf. Figure[2j). A g(b,a) in the (6, a)-plane color-coded as 
specified in color bar for all three taxation schemes. Solid lines divide zones of income growth from 
income destruction. Dashed lines are optimal tax rates for given admin rate a op t(b). Above the dotted 
line income destruction must happen. B Lines of A in one plot for comparison. Colors indicate tax 
schemes. Linestyles as above. 
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Figure 4. Optimal tax rate, maximal growth rate and maximal income of government. A 

Optimal tax rate a op t{b), B maximum achievable growth rate g m ax(&), C governmental income rate 
under optimal tax rate gov opttax (&) for different tax schemes: proportional tax (blue), regressive tax 
(green), progressive tax (red). Stars indicate the location of the maximum of gov opttax in all plots. The 
dashed line is to show that under the maximal gov opttax proportional taxation gives the largest growth 
factor. Parameters as in[3l 
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Figure 5. More agents, less and more risky human capital. Extension of Figure |3jA. with less 
and more risky human capital production functions (in rows) and more agents (in another column). 




Figure 6. Zones of growth and destruction for more agents, less and more risky human 
capital. Simulation results analog to Figure [3j3 with less and more risky human capital productions 
functions (in rows) and more agents (in another column). (Some lines as in Figure P7[) 
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Figure 7. Zones of growth for different population sizes. Same lines as in Figure |6j but such 
that different population sizes are in one plot and taxation schemes in columns. 
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Figure 8. More agents, less and more risky human capital and its impact on optimal tax 
rate, maximal growth rate and maximal income of government. Simulation results analog to 
Figure HI with less and more risky human capital production functions and more agents. 



Supporting Material 

matlab-code to run simulation and produce figures: 

function figs 

'/, Produces all simulation data and all figures for the Paper 
'/, by Jan Lorenz, Fabian Paetzel and Frank Schweitzer, 2012 
'/. programmed by Jan Lorenz 

Fig 1; 4 subfigures: figExRedis f igExRedisdynf ee f igExRedisproptax f igExRedisdynmax 
exampleRedistribution 

% Fig 2; 4 subfigures: f igMSPgrowthExMain f igMSPgrowthExSubl f igMSPgrowtnExSub2 f igHSPgrowthExSub3 
exampleTrajectories 

'/, Generate systematic simulation data 

dataHSPgrowthrates 7. WARNING TAKES TIME! PROBABLY AT LEAST A WEEK ... 
"/. Figs 3—8; 

colmap '/, save the colormap skewjet to file 
figColormap figure of colormap 

f igSimulation calls subfunctions plotf igs_f ixedN, plotf igs_taxf unction and produces lots of subfigures 



FUNCTIONS PRODUCING FIGURES 
function exampleRedistribution 
figure{l);clf 

y = [100 300 600 1000 1500 2100] ' ; 

a = 1/3; 

b = 0.25; 

f = getfee(y,a) ; 

m = getmax(y,a) ; 

ylim([0 l.l*max(y)]);xlim([0 length(y)+l] ) ;hold on 

xfill=[0.7 0.7 1.3 1.7 2.3 2.7 3.3 3.7 4.3 4.7 5.3 5.7 6.3 6.3]; 

yfill=a*[0 yCl) y(l) yC2) y(2) y(3) y(3) y(4) y(4) y(5) yC5) y(6) y(6) 0]; 

f ill{xf ill,yfill, [0.5 0.5 1] , 'EdgeColor' , [0 1] , 'LineWidth' ,2, 'FaceAlpha' ,0.5) ; 

fill([0 length(y)+l length(y)+l 0] , [m m max(ylim) max(ylim)] , [1 0.5 . 5] , ' LineStyle ' , ' none 1 , 'FaceAlpha ' , . 5) ; 

fill([0 length(y)+l length(y)+l 0] , [f f 0] , [0 . 5 1 . 5] , 'LineStyle ', 'none ',' FaceAlpha' ,0 . 5) ; 

plotf [0 length Cy)+1] , [m m] , 'r ' , 'Linewidth' ,2) ; 

plot([0 length(y)+l] , [f f] , 'g" , 'Linewidth' ,2) ; 

bar(y, 'FaceColor' ,0.25*[1 1 1] , 'BarWidth' ,0.6) ;hold on; 

set(gca, 'Layer* , 'top' , 'Xtick' ,l:length(y) , 'Position' , [ 0. 1600 0. 1500 0.7750 0.8150]) 
make_f iguresC 'figExRedis ' , [82 192 215 332] , 'removeEps ' , 1) 
mode = -['proptax' , 'dynfee' , 'dynmax'}; 
for k=l:3 



afterRedis(k+l,mode{k},y,a,b) 

make_f igures ( ['f igExRedis ' mode{k>] , [5 243* (k-1) . 6* [269 221] ] , 'removeEps ' ,1) 



function exampleTraj ectories 
tmax=500; 

7, parameters for 4 figures 
n = [10 10 10 100] ; 
tax = [0.3 0.3 0.01 0.3] ; 
adm = [0.2 0.6 0.2 . 2] ; 

log-normal 
Mean = 1.5; geoHean = 1/Mean; 
mu = logCgeoHean) ; si = sqrt(2*(log(Mean)-mi 
r = exp(mu + si*randn(n(end) ,tmax) ) ; 
maxrate=Mean* ( 1-adm . *tax) ; 
for i=l:length(n) 
figure(i);clf 

x = ones(n(i) ,tmax) ;xcol='g' ; 
y = ones(n(i) ,tmax) ;ycol='b' ; 
z = ones(n(i) ,tmax) ;zcol='r' ; 
for t = 2:tmax 



x(:,t) = 

y(=,t) = 

z(:,t) = 

plot(sum(x)' 
plot (sum(y) ' 
plot(sum(z)' 



redis(prod(HCprod(x( 
redis(prod(HCprod(y( 
redis(prod(HCprod(z( 

:ol) ;set(gca, 'YScale' 
ol) ;set(gca, 'YScale' 
ol) ;set(gca, 'YScale' 



,t-l), 
,t-l), 
,t-l), 



randomvecto: 
randomvecto: 
randomvecto: 



n(i) ,t) ) ) , 'dynf ee' , 'tax' ,tax(i) , 'adm' , adm(i) ) ; 
n(i) ,t))) , 'proptax ' , 'tax' ,tax(i) , ' adm' , adm(i) ) ; 
n(i) ,t))) , ' dynmax ' , ' tax' ,tax(i) , ' adm' , adm(i) ) ; 



'log' ) ;hold < 

'log') 

'log') 



plot (n(i)*maxrate(i) . " (1 :tmax) , 'k' ) ; 
pi = n(i)*geoMean. "(iitmax) ;plot(pl(pl>0) , 'k') ; 



set(gca, 'Posit: 
text ( 1 . 01 *tmax , sum (x ( 
text ( 1 . 01 *tmax , sum (y ( 
text (1.01*tmax, sum (z( 
text (1.01*tmax, (n(i) *maxrate (i) ~tmax) 
xlabel('t') 

ylabel (' total income (\Sigma_i y_i(t) ) ' ) ; 



[0.14 0.1 0.65 0.8]) 

end) ) , ' regr . tax' , ' color ' ,xcol) 
end) ) , 'prop, tax' , 'color' ,ycol) 
end)),'prog. tax' , 'color' ,zcol) 

N((l-ba)<\eta>)~t' 



'color' , 'black' ) 



switch : 



title(['N = ' num2str(n(i) ) ', a = ' num2str(tax(i)) ', b = ' num2str (adm(i) )] ) ; 
case 2 

title(['N = ' mim2str(n(i)) ', a = ' num2str(tax(i)) ', {\bf b = ' num2str(adm(i)) '}']); 
case 3 

title(['W = ' num2str(n(i)) ', {\bf a = ' num2str(tax(i) ) '>, b = ' num2str (adm(i) )] ) ; 
case 4 

title(['{\bf N = ' num2str(n(i)) '}, a = ' num2str(tax(i) ) ', b = ' num2str (adm(i) )] ) ; 
xlim([0 500]) 

ylim( [min( [sum(x) sum(y) sum(z)] )*10~-10 max(n(i) *maxrate (i) "tmax) ~1 . 05] ) 



g=axes('position',[0.19 0.65 0.17 0.15]); 
X = 0:0.01:5;plot(X,lognpdf(X,mu,si),'k'); 

set(g, 'XTick' , [0 1 2] , 'YTick' , [] , 'YColor' , 'w' , 'Clipping' , 'off ') ; 
ylim([0 .9*max(ylim)] ) 

text (0.8,0. 9*max(ylim) , [ ' \eta log-normal ' ] ) 

text(1.2,0.6*max(ylim) , ['\mu=' num2str(round(mu*1000)/1000) ' , \sigma=' num2str(round(si*1000)/1000)] ) 

text(1.7,0.3*max(ylim) , ['<\eta> = ' num2str(round(Mean*1000)/1000) ', <\eta>_{geo> = ' num2str (round (geoHean*1000) /1000)] ) 
title('HCprod, pdf \eta') 



switch : 



make_f iguresf 'f igMSPgrowthExHain' 
case 2 

make_f igures('f igMSPgrowthExSubl ' 
case 3 

make_f igures ( 'f igMSPgrowthExSub2 ' 
case 4 

make_f igures ( 'f igMSPgrowthExSub3 ' 



[877 373 400 350] ) ; 



[600 
[800 
[1000 



[400 
[400 
[400 



350] /2.3]); 
350J/2.3]); 
350] /2.3]); 



function afterRedis(k,mode,y,a,b) 
figure(k);clf ; 
pg=sum(y)*a*(l-b) ; 
switch mode 

case 'proptax' ; lcol = [0.5 0.5 l];fcol =[0 1]; 
case 'dynfee' ;lcol = [0.5 1 0.5] ;f col =[010]; 
case 'dynmax' ;lcol = [1 0.5 0.5] ;f col = [10 0]; 

bar (1 : length (y) ,redis (y ,mode , ' adm ' ,b , ' tax ' , a) , 'FaceColor ' ,f col , 'BarWidth ' ,0.5) ; hold on 
bar(l : length (y) ,ones (1 , length (y) ) *pg/length(y) , 'FaceColor' ,lcol, 'BarWidth' ,0.5) ; 
set (gca, 'Layer',' top', 'Xtick ' , 1 : length(y) , ' Position', [ 0.200 0.1500 0.7750 0.8150]) 
box on;ylim([0 1 . l*max(y)] ) ;xlim( [0 length(y) +1] ) ; 



function figColormap 
figure(l);clf 
load colmap 

c = exp([-1.5 l]);caxis(c); 

imagesc(linspace(c(l),c(2)) , 1 , linspace (c (1) ,c(2))); 
colormap ( skewj et ) 

set(gca, 'YTick' ,[], 'XAxisLocation' , 'top' , 'Position' , [0. 1 0.35 0.8 0.15]) 
text (0.99,2.5, ' | growth \rightarrow ' ) 

text (0.99,2.5, ' \left arrow destruction ' , 'HorizontalAlignment ' , 'right ' ) 
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title ({' color axis: average growth factor g{b,a)0) 
make_f igures( 'f igColormap' , [508 500 313 100] ) 

function f igSimulation 

mean_eta = [1.25, 1.5, 3]; geomean.eta = l./mean_eta; 
MU = log(geomean_eta) ; SI = sqrt(2*(log(mean_eta)-MU)) ; 
for n = [10 100] 

for dist.ind =[12 3] 

mu = round(100*MU(dist_ind)) ; si = round(100*SI (dist.ind) ) ; 

plotf igs_fixedN(n,mu, si) 

modes = {'dynfee 1 , 'proptax' , 'dynmax' }; 
for mode = 1:3 

for dist.ind =[12 3] 

mu = round(100*MU(dist_ind)) ; si = round(100*SI (dist_ind) ) ; 

m = [10 100] ; 

plotf igs_taxf unction (modes{mode},mu, si, NN) 

function plotf igs_taxf unction (taxf unction, muName,siName,NN) 
modes = [] ; mu= [] ; '/, necessary initialization 
figure{l);clf ; 
for nlnd = l:length(NN) 

load(['data' 'N' num2str(NN(nInd) ) 'mu' num2str(muName) 'si' num2str(siName) '.mat']) 

mode = find(strcmp(taxf unction , modes) ) ; 

if length(NN) == 2 

linecol = ones (2 ,3) - [0 . 33* [1 11]; [111]]; 

linecol = ones (3,3) - [0 . 33* [1 1 1]; 0.66*[1 11]; [111]]; 

switch mode 

case 1; linecol(: ,2)=1; 
case 2; linecol(: ,3)=1; 
case 3; linecol (:, 1) =1 ; 

g=growrate(l,l, 1, : ) ;g=g( : ) ;growrate=growrate ( : , : , : , ~isnan(g) ) ;f ind(~isnan(g) , 1 , 'last') 
P = mean(growrate( : , : ,mode , : ) ,4) ; 
[tmp,maxind]=max(P) ; argmaxA = A(maxind) ; 

[C(nlnd) .h.h(nlnd) .h]=contour(B, A,P, ' color' .linecol (nlnd, : ) , 'linewidth' ,2, 'LevelList' ,0) ; hold on 
T(nlnd) .T=clabel(C(nInd) .h.h(nlnd) .h, 'LabelSpacing' ,500, 'BackgroundColor ' , [1 1 1] ) ; 
set (T(nlnd) ,T, 'string' ,num2str(NN(nInd) ) ) ; 
plot (B , smooth ( argmaxA) , ' — ' , ' color ' , [linecol (nlnd, : )] ) ; 

f igure(l) ; 

axis equal;axis( [min(B) max(B) min(A) max(A)]) 
switch mode 

case 1; titlestr = 'regressive tax'; 

case 2; titlestr = 'proportional tax'; 

case 3; titlestr = 'progressive tax'; 

plot(B,l./B*(l-l/exp(mu+si-2/2)) , 'k: ') ; 

make_f igures( ['f igSim' taxfunction 'mu' num2str (muName) 'si' num2str(siHame)] , [360 530 233 193]) 
function plotf igs_fixedN(n,mu, si) 

load( [ ' data' ' N ' num2str (n) 'mu ' num2str (mu) ' si ' num2str (si) ' . mat ' ] ) 

g=growrate( 1,1,1, : ) ;g=g( : ) ; growrate=growrate ( : , : , : , ~isnan(g) ) ;f ind(~isnan(g) , 1 , 'last') 

f igure(l) ; elf ; figure (2) ; elf ; figure(3) ; elf ; 

linespec={ 'g' , 'b' , 'r'>; 

maxadminlncomelnd = zeros (1,3); 

maxmaxP = zeros(l,3); 

for mode=l:3 

P = exp(mean(growrate ( : , : .mode , : ) , 4) ) ; '/, mean(growrate ( : , : .mode , : ) ,4) is log of growthrate of model y (t) =N*exp( . . . ) ~t 
% P is thus geometric mean of growthrates y(t) = N*P"t 
[maxP ,maxind] =max(P) ; argmaxA = A(maxind); 
smoothargmaxA=smooth(argmaxA) ' ; 
adminlncome = smooth (maxP . *smoothargmaxA . *B) ' ; 
[tmp , maxadminlncomelnd (mode) ] = max (adminlncome) ; 
maxmaxP (mode) = maxP(maxadminIncomeInd(mode)) ; 
figure (1) 
subplot ( 1 ,3 ,mode) 

load( ' colmap 1 ) ; '/, load the colormap skewjet, can be replaced by a default one 
colormap ( skewj et ) 

imagesc (B , A,P) ; axis xy;axis equal;axis tight; 
hold on; 

contour(B,A,P, [1 1] ,' color ', 'k ' ) 
plot (B , smoothargmaxA , 'k — ' ) 
c = exp([-1.5 1] ) ;caxis(c) ; 

if mode==l;ylabel('tax rate a');end; if mode==2 ; xlabel (' admin rate b');end 
switch mode 

case 1; title('{\bf regressive tax}' , 'color' , 'g') ; 
case 2; title('{\bf proportional tax}' , 'color' , 'b') ; 
case 3; title('{\bf progressive tax} ', ' color ', 'r' ) ; 

plot(B,l./B*(l-l/exp(mu+si~2/2)) , 'k: ') ; 
figure (2) ; 

contour(B,A,P, [1 1] ,' color ', linespec{mode} ,' linewidth' , 2) ; hold on 
C(mode)=plot( [100] , [100] , linespec-frnode} , 'linewidth' ,2) ; '/, just for the legend 
plot (B , smoothargmaxA, [linespec{mode} ' — ' ] ) ; 
figure (3) ; 
subplot (3, 1,1) 

plot (B , smoothargmaxA, [linespec{mode} ' — ' ] ) ; hold on 

plot (B (maxadminlncomelnd (mode) ) , smoothargmaxA (maxadminlncomelnd (mode) ) , [linespec{mode} ' * '] ) 
subplot(3,l,2) 

plot([0 0.8],[1 1] , 'color' ,0.5*[1 1 l]);hold on 
plot (B ,maxP , 'color' , linespec{mode}) ; 

plot (B (maxadminlncomelnd (mode) ) ,maxP (maxadminlncomelnd (mode) ) , [linespec{mode} '*']) 
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subplot (3, 1,3) 

plot (B .adminlncome , 'color ' , linespec{mode}) ; hold on 

plot (B (maxadminlncomelnd (mode) ) , adm in Income (maxadminlncomelnd(mode) ) , [line spec {mode} '*']) 
f igure(l) 

make_f igures( ['fig' 'Sim' ' 1 ' 'W num2str(n) 'mu' num2str (round (mu* 100) ) 'si' num2str (round (si* 100))] , [6 533 522 188] ) 
f igure(2) ; 

axis equal;axis( [min(B) max(B) min(A) max(A)]) 
plot(B,l./B*(l-l/exp(mu+si~2/2)) , 'k: ') ; 
ylabelf'tax rate a' ) ; xlabel ( ' admin rate b'); 

make_f igures( ['fig' 'Sim' '2' 'N ' num2str(n) 'mu' num2str (round (mu* 100) ) 'si' num2str (round (si* 100))] , [985 435 292 288] ) 
f igure(3) 

subplot (3,1,1) ; '/.title ({'optimal tax rate'}) 
ylabel('a_{opt}');xlim([0 . 8] ) ;ylim ( [0 1]) 
subplot(3,l,2) ; '/.title ( 'maximal growth factor') 

ylabel('g_{max}');axis tight ;yax=ylim; xlim( [0 . 8] ) ; ylim( [0 . 9*yax(l) 1 . 05*yax(2) ] ) 
[tmp.i] = max(maxmaxP) ; 

plot([0 0.8] ,max(maxmaxP)*[l 1] , [linespec{i> ' — ']) 
subplot(3,l,3); 7,title('government income') 

xlabeK'admin rate b' ) ;ylabel ( 'gov_{opttax> ' ) ; axis tight ;yax=ylim ;xlim( [0 . 8] ) ; ylim( [0 1 . 05*yax(2)] ) 

make_f igures( ['fig' 'Sim' '3' 'N ' num2str(n) 'mu' num2str (round (mu* 100) ) 'si' num2str (round (si* 100))] , [808 156 213 634] ) 



7, FUNCTIONS PRODUCING DATA 
function dataMSPgrowthrates 

% makes samplesize=100 simulation runs for each data point for different 

7. values of mu, sigma, W, and taxregime; and estimate a growthrate for 

7 t each; writes this data to files to be read be figure-producing functions 

mean_eta = [1.25, 1.5, 3]; geomean_eta = l./mean_eta; 

HU = log(geomean_eta) ; SI = sqrt (2* (log (mean_eta) -MU) ) ; 

N = [10 100] ; 

for n = N 

for dist_ind = l:length(mean_eta) 

mu = MU(dist_ind) ; si = SI (dist_ind) ; 

files = dir; 

if ~any(strcmp({f iles .name}, ['data' 'N' num2str(n) 'mu' num2str(round(mu*100)) 'si' num2str(round(si*100) ) '.mat'])) 
7, external parameter 
tmax=500 ; 
samplesize = 100; 
7, running parameters 
A = 0:0.02:1; 
B = 0:0.02:0.8; 

modes = { 'dynf ee ' , 'proptax ' , ' dynmax 1 > ; 
y = ones(n,tmax) ; 

growrate = nan(length(A) ,length(B) ,length(modes) , samplesize) ; 
SSE = nan(length(A) ,length(B) .length(modes) , samplesize) ; 
RMSE = nan(length(A) ,length(B) .length(modes) , samplesize) ; 
RSqUARE = nan(length(A) ,length(B) , length(modes) .samplesize) ; 
else '/, when a file exists one might resume and finish an old computation 
tmax = 500; 
samplesize = 100; 
y = ones(n,tmax) ; 

load(['data' 'N' num2str(n) 'mu' num2str (round (mu*100) ) 'si' num2str(round(si*100)) '.mat']) 

for s=l : samplesize 

r = exp(mu + si*randn(n , tmax) ) ; 
for a=l:length(A) 

for b=l:length(B) 

for m=l : length (modes) 

if isnan(growrate(a,b,m,s)) 

y(:,D = i; 

for t = 2: tmax 

y( : , t) = redis(prod(HCprod(y( : ,t-l) , 'randomvector ' ,r ( : ,t))) ,modes{m>, 'tax' ,A(a) , 'adm' ,B(b)) ; 

'/, fit for polyl: ly = pl*ly + p2 

ly = log(sumCy))'; 

if all(~isinf (abs(ly))) 

[f.gof] = fit((0:tmax-l) ' ,ly, 'polyl' , 'Lower' , [-inf log(n)] , 'Upper' , [inf log(n)]); 

growrate(a,b,m,s) = f.pi; 

SSE(a,b,m,s) = gof.sse; 

RMSE(a,b,m,s) = gof.rmse; 

RSQUARE(a,b,m,s) = gof.rsquare; 

save( ['data' 'N ' num2str (n) 'mu' num2str (round (mu* 100) ) 'si' num2str (round (si*100)) ' .mat'] , 'A' , 'B' , 'RMSE' , ' RSQUARE ' , 'SSE' , 'growrate' , 'modes' , 'mu' , 'si' , 'tmax') 



7. FUNCTIONS FOR COMPUTATION 



function y = prod(h, varargin) 

y - hi 

function h = HCprod(y, varargin) 

% y is a row vector of the current wealth of n=length(y) agents 
7. defaults and varargins 
r - [] ; 

for i=l : length (varargin) /2 
switch varargin{i*2-l} 

case 'randomvector' ;r = varargin{i*2}; 



7, production mechanism 
h = y.*r; 



function yn = redis(y, mode, varargin) 
7, defaults and varargins 
tax = [] ; 
adm = [] ; 

for i=l : length (varargin) /2 
switch varargin{i*2-l} 

case 'tax';tax = varargin{i*2]*; 
case 'adm';adm = varargin{i*2]*; 

7 t redistribution mechanism 
switch mode 

case 'proptax' 

7, public good, prop tax, admin cost 

if isempty(tax) ; tax = 0.5; end "/, tax rate for the public good 
if isempty(adm) ; adm = 0; end 7« rate of administrative cost for public good 
t = sum(y*tax) ; 
pg = sum(t)*(l-adm) ; 
yn = y*(l-tax)+pg/length(y) ; 
case 'dynfee' 

% public good, a fixed amount is taxed (or all) 

if isempty(tax) ; tax = 0.5; end 7« tax rate, what proportion to get from total wealth 
if isempty(adm) ; adm = 0; end 7« rate of administrative cost for public good 
f = getfee(y,tax) ; 
t = min( [y f *ones (size(y) ) ] , [] ,2) ; 
pg = sum(t)*(l-adm) ; 
yn = (y-t)+pg/length(y) ; 
case ' dynmax' 

X public good, all over max is gathered for pg, admin cost 

if isempty(tax) ; tax = 0.5; end 7« tax rate, what proportion to get from total wealth 
if isempty(adm) ; adm = 0; end % rate of administrative cost for public good 
m = getmax(y,tax) ; 

t = max( [zeros (size (y) ) y-m],[],2); 

pg = sum(t)*(l-adm) ; 

yn = (y-t)+pg/length(y) ; 



function f = getf ee(y,tax) 

X y is row vector of numbers, tax is a desired tax rate 

7, extracts f such that sum (max( [zeros (size (y) y-m] , [] ,2))= sum(y)*tax 

y = [0 ; sort(y, 'ascend')] ; 

t = sum(y)*tax; 

n = length(y) ; 

i = i; 

while sum(y(l:i-l))+y(i)*(n-i+l)<t 
i=i+l; 

d = sum(y(l:i-l))+y(i)*(n-i+l)-t; 
f = y(i) - d/(n-i+l); 

function m = getmax(y ,tax) 

X y is row vector of numbers, tax is a desired tax rate 

7 t extracts m such that sum(max( [zeros (size (y) y-m] , [] ,2))= sum(y)*tax 

y = [sort(y, 'descend') ; 0]; 

t = sum(y)*tax; 

i = 1; 

while sum(y (1 : i)-y (i+1) ) <t % optimizing potential for large y! 
i=i+l; 

d = sum(y(l:i)-y(i+l))-t; 
m = y(i+l) + d/i; 

X MISCELLANEOUS FUNCTIONS FOR PRINTING AND C0L0RMAPS 

function make_f igures (savestr ,pos .varargin) 
% defaults removeEps = 0, pdfcrop = 
removeEps = 0; 
pdfcrop = 0; 

if mod (length (varargin) ,2)==1 

error (' varargin not correct'); 

else 

for i=l : length (varargin) /2 
switch varargin{i} 
case 'removeEps ' 

removeEps = varargin-[i+l]- ; 
case 'pdfcrop' 

pdfcrop = varargin{i+l]*; 

set (gcf , ' PaperPositionMode ' , ' auto ' , ' Position' , pos) 

print (savestr , ' -deps2c ' ) 

eval( [ ' ! epstopdf ' savestr ' . eps ' ] ) 

if removeEps==l 

eval ( [ ' ! rm ' savestr ' . eps ' ] ) 

if pdfcrop==l 

eval ( [ ' ! pdfcrop ' savestr ' .pdf ' savestr ' . pdf ' ] ) 



function colmap 

7 t Code ommitted. Just saves a particular color map as a file to load it for graphics 



